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Abstract 

This  paper  deals  with  the  diffraction  of  a  sound 
pulse  due  to  an  instantaneous  line  source  by  a  parallel, 
rigid,  fixed  circular  cylinder.   This  problem  can  also 
be  interpreted  in  electromagnetic  terms.  The  diffracted 
field  is  resolved  into  terms  representing  diffracted 
pulses  which  have  encircled  the  cylinder  a  niiinber  of 
times  in  either  sense,  and  these  in  turn  are  expressed 
as  a  series  of  'propagation  modes'  with  the  same  pulse 
fronts.  Approximations  valid  near  the  front  of  a  mode 
are  then  obtained.  The  corresponding  results  for  an  in- 
cident plane  pulse  are  derived  by  a  limitir^  process.   It 
is  found  that  a  diffracted  pulse  has  an  essential  singu- 
larity of  a  definite  type  et  its  fronts,  similar  to  the 
initial  stages  of  diffusion,  and  that  the  pressure  rises 
very  gradually  et  first.  But  after  time  intervals  of  the 
order  of  one-third  to  one-half  of  the  radius  of  the  cylin- 
der divided  by  the  velocity  of  sound,  there  seems  to  be 
no  appreciable  diffraction  effect. 
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Introductlon 

The  pressvire  field  behind  the  front  of  a  sotmd  pulse  can  be  discussed 
quite  easily  if  it  is  assumed  that  it  is  a  regular  function  of  the  time  elapsed 
since  the  arrival  of  the  piolse.  The  front  itself  is  propagated  according  to  the 
laws  of  geometrical  optics j  the  square  of  the  pressure  jump  at  the  front  obeys 
the  intensity  law  of  geometrical  optics  *-  *  -"j  finally,  the  derivatives  of  the 
pressiire  with  respect  to  time  at  the  front  can  be  calculated  recursively,  so  that 
the  pressiore  can  be  e3q)anded  in  an  ascending  series  of  powers  of  the  time  counted 
from  the  arrival  of  the  front. 

When  a  sound  pulse  is  scattered  by  an  obstacle,  this  method  can  be  ap- 
plied behind  a  reflected  front,  but  not  behind  a  diffracted  one.  To  distinguish 
between  these  two  types  of  front,  consider  for  simplicity  a  convex  obstacle  which 
is  fixed  and  rigid,  and  an  incident  pulse  with  an  expanding  front.  The  origin  of 
time  can  be  chosen  at  an  instant  when  the  obstacle  is  in  the  region  not  yet  dis- 
turbed by  the  incident  pulse.  The  surface  of  the  obstacle  can  be  divided  into 
two  domains:  the  reflecting  region,  every  point  of  which  can  be  joined  to  the 
initial  position  of  the  incident  front  by  an  incident  ray  in  the  exterior  of  the 
obstacles,  and  the  shadow.  The  reflected  fronts  are  the  surfaces  normal  to  the 
reflected  rays,  Wiich  issue  from  the  reflecting  region  on  the  obstacle  according 
to  the  reflection  law  of  geometrical  optics.  In  the  shadow,  the  disturbance  pro- 
pagates over  the  obstacle  with  the  speed  of  sound  and  is  generally  bounded  by  dif- 

[2] 
fracted  fronts  which  are  involutes  of  the  surface  of  the  obstacle  ■-  -• .  The  ex- 
pansion of  the  pressure  in  ascending  powers  of  the  time  counted  from  the  arrival 
of  the  front  fails  at  a  diffracted  front  because  the  surface  of  the  obstacle. 
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where  the  boundary  conditions  must  be  applied,  is  a  caustic  of  the  system  of  dif- 
fracted fronts. 

It  is  therefore  of  some  interest  to  investigate  the  behavior  of  the 
pressure  behind  a  diffracted  front  in  a  specific  case*  The  diffraction  problems 
with  knovm  solutions  -  such  as  the  half-plane  and  wedge  problems  -  are  instances 
where  diffraction  is  essentially  an  edge   effect.  The  simplest  problem  involving 
a  convex  obstacle  is  that  of  the  diffraction  by  a  circular  cylinder  of  two-dimen- 
sional incident  pulses  in  which  the  pressure  is  the  same  at  all  points  of  any 
straight  line  parallel  to  the  axis  of  the  cylinder.  This  is  the  problem  dis- 
cussed in  this  paper,  attention  being  concentrated  on  the  diffracted  pressure 
field*. 

The  problem  which  is  here  formulated  and  solved  in  terms  of  sound  pulses 
can  also  be  interpreted  in  electromagnetic  terms.  Suppose  that  the  cylinder  is  a 
perfect  conductor,  and  that  the  incident  field  is  such  that  the  only  non^^anishing 
component  of  the  magnetic  field  is  that  parallel  to  the  axis  of  the  cylinder,  '^hen 
it  satisfies  the  same  boundary  conditions,  (1.6),  at  the  cylinder,  as  the  pressure 
does  in  the  acoustic  case.  The  other  two-dimensional  electromagnetic  problem,  in 
which  the  electric  field  is  parallel  to  the  axis  of  the  cylinder,  can  be  solved  by 
the  same  technique;  the  modifications  required  are  trival,  and  some  of  the  calcula- 
tions are  actually  simplified. 

1.  The  Green's  Function 

To  simplify  the  work,  the  radius  of  the  cylinder  will  be  taken  as  the 
unit  of  length,  and  the  velocity  of  sound  as  the  unit  of  velocity.  Then. the  unit 


A  related  problem  is  that  of  the  propagation  of  a  sound  pulse  in  a  stratified 
atmosphere  bounded  by  a  rigid  plane,  in  vhich  the  sound  velocity  decreases  with 
distance  from  the  bounding  plane.  The  case  in  which  the  reciprocal  of  the  sound 
velocity  increases  linearly  with  the  height  has  been  investigated,  and  an  account 
of  it  will,  it  is  hoped,  be  published  as  a  separate  paper |  by  the  author. 
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of  time  is  the  time  taken  hj  a  pulse  front  to  travel  a  distance  equal  to  the 
radius.  Only  two-dimensional  problems  viill  be  considered.  With  the  axis  of 
the  cylinder  as  the  z  axisj  the  state  of  the  medium  is  completely  specified 
by  giving  the  excess  pressvire  p  and  the  particle  velocity  components  u,  v,  as 
functions  of  x,  y  and  the  time  t.  The  equations  of  sound  propagation  are 

where  K  is  the  acoijstic  impedance  (the  product  of  density  and  sound  velocity)  of 
the  medium,  and  a  ••  source"  term  q  has  been  included  in  the  equation  of  continuity 
(1.2)  as  this  is  sometimes  a  convenient  way  of  giving  the  incident  pulse.  •^  velo- 
city potential  0  can  now  be  defined  by  the  equations 

(1.3)  p  ■  i.  >=''  =  -i.  fr  -  -^  • 

It  satisfies  the  wave  equation 

(l.Ii)  1|    -     v2  jzf     =     q     . 

at 

The  pressure  p  also  satisfies  the  wave  equation,  with  a  forcing  term  g^,  and  one 
can  work  with  either  0  or  p. 

If  (r,6)  denote  plane  polar  coordinates  defined  as  i;isual  by 

(1.5)  X  =  r  cos  e,    y  =  r  sin  0 

then  the  boundary  condition  at  the  cylinder,  which  expresses  the  vanishing  of  the 
normal  component  of  the  particle  velocity,  is 


(1.6)  [221 


dr  Jr  =  1 


-U- 


</-oo 


r  =  1 


An  incident  pulse  may  be  specified  in  variovis  ways.  It  may  be  given 
explicitly;  if  the  origin  of  time  is  taken  at  an  instant  prior  to  that  when  the 
incident  fraitto\iches  the  cylinder  for  the  first  time,  this  is  equivalent  to 
giving  0  and  -^    when  t  ■  0»  Alternatively,  these  initial  values  nay  be  given; 
they  are  equivalent  to  a  distribution  of  instantaneous  sources  and  doublets  on 
the  plane  t  «  0  in  space  time.  Finally,  a  source  distribution  q  which  generates 
the  incident  pulse  may  be  specified  directly.  In  all  cases,  the  diffraction 
problem  may  be  solved  if  its  Green's  function  is  known.  This  may  be  defined  as 
the  solution  of 

2 
(1.7)  L|.  -  V^G  =  2n  5(x  -  x  )6(y  -  y  )5(t  -  t  ) 

dt 

which  satisfies  the  boimdary  condition 

(1-8)  [i]r<    -    ° 

(here,  as  elsewhere,  5  denotes  the  Dirac  5-function) .  The  general  case  can  ob- 
viously be  deduced  by  rotation  and  translation  from  the  special  case  x  =  r  , 

y  =  t  =  0,  so  that 
•'o    o    ' 

2 

(1.9)  L^  .  T^G  -  2n  5(r  -  r  )5(y)5(t) 

at  ° 

and  this  alone  will  be  discussed  here.  If  the  cylinder  were  absent,  the  solution 
of  (l»9)  would  be  Hadamard's  elementary  solution 

(1.10)  G^(x,y,z}r^)  -  (t^  -  R^)"^/^  H(t  -  R) 
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where  H  is  Heaviside's  unit  function  and 

(1.11)  R  -  J(x  -  r^)2  *  y^T'^^ 


This  is  the  incident  field  which  generates  the  Green's  function.  (It  must  be 
understood  throughout  that  r  >  1) • 

It  is  natural  to  introduce  the  polar  coordinates  (1.5).  One  can  then 
think  of  G  as  defined  not  only  in  the  "  physical  plane"  (r  >  0,  -n  <  6  <  n)  but 
also  on  a  Rieraann  surface  whose  sheets  su:e  given  by  (2m  -  l)n  <  6  <  (2m  +  l)n, 
m  =  ...j  -1,  0,  1,  ...  •  Obviously,  the  instantaneous  soiirce  atx=>r,y  =  0 
then  corresponds  to  an  infinity  of  sources  situated  at  the  points  r  ■  r  ,  0  =  2nin 
(m  =  ••.  -1,  0,  1,  ...  ).  Thus  (1.9)  mxist  be  replaced  by 

(1.12)  ^  -  V^G  =  |i  6(r  -  r^)5(t)  21   5(e  +  2m)    . 
dt  o  m=-  00 

(The  factor  1/r  is  required  because  the  element  of  area  of  polar  coordinates  is 
rdrd©).  If  follows  at  once  that  one  may  put 

00 

(1.13)  G(r,©,tir  )  »  2Z   ^(^>®  *  2mn,  t^  r  ) 


m=-oo 


where  F  is  defined  for  -co  <  6  <  oo  by  the  equation 


2 

(I.IU)  ^  -  V^F  -  ^  6(r  -  r^)5<@)§<fe) 

dt  o 


and  the  boundary  condition 
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One  may  think  of  F  as  the  field  due  to  an  incident  pulse  generated  by  a  modified 
"source"  which  has  the  requisite  singularity  only  in  the  sheet  |©|  <  n  of  the 
Riemann  surface.  The  formula  (1.13)  is  thenan  instance  of  the  general  reciprocity 
theorem. 

2.  The  Laplace  Transform  of  the  Green's  Function 
Let 

CD 

(2.1)  F(r,e,S}r^)   =   /  F(r,e,t)e"^*  dt 

be  the  Laplace  transform  of  F,  and 

-00 


(2.2)         F*(r,v,sjr  )  =  ^W-   /  F(r,e,s}r^)e"^^®de 

°      (2n)^/2  /_^         ° 


be  the  Fourier  transform  of  F.  Then*" 

(2.3)  s^F  -  V^F  =  —  6(r  -  r  )5(e) 

o 

and 

Also,  it  follows  from  (l.l5)  that 


^ 

■^  Since  it  is  known  that  F  is  singmLar  on  certain  characteristxcs  in  space-time,  the 
statement  that  (2.3)  is  the  Laplace  transform  of  (1.13)  is  open  to  doubt.  One  can 
prove  its  validity  by  using  the  concept  of  the  "  finite  part"  of  a  divergent  inte- 
gral due  to  Hadamard,  or  by  interpreting  F  as  a  "  distribution"  in  the  sense  of 
L.  Schwartz. 
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It  may  be  assumed  that  F  is  continuous  at  r  ■  r  .  Then  (2«ii)  is  equivalent  to 


{2.6) 

and 

(2.7) 


d^F*  ^  1  dF*   /v^    2\  *     r.       I      J       ^ 


r  +  0 
0 

r  -  0 

L  J  O 


0, 


dF 
dr 


r  +  0 

0 

r  -  0 
J  0 


(2n) 


1/2 


Adding  the  condition  F  -^0  as  r  ->oo,  one  finds  after  a  short  calculation  that 

^       K.(sr  ) 


(2.8)  F  (r,v,sjr^)  -  I 


,(sr)    r  ,  ,  n 

T \l   (sr  )K  (s)  -  I  (s)K  (sr  )l  ,  (r  >  r  )  , 

(s)   I  '^   O  V       V    V   O    {  *  -  o   ' 


1/2  \^^^^ 
K 


where  I  ,  K  denote  the  modified  Bessel  functions  of  order  v  of  the  first  and 
second  kind  respectively.  That  (2,8)  satisfies  (2.$)  is  obvious;  the  validity 
of  (2.7)  follows  from  the  identity 


(2.9) 


i;(z)K^(z)  -  i,(z)k;(z)  -  \ 


The  Laplace  transform  F  can  now  be  obtained  at  once  in  the  form 


(2.10)  F(r,e,sjr  )  = 

'    '  o 


'^     I  (sr)K  (s)  -  I  (s)K  (sr) 

V    V      V    V     „  /   >,  ive  , 

K  (sr  )e    dv   , 


-00 


Vs) 


V    0 


for  r<r:ifr>r,  then  r  and  r  must  be  interchanged. 
—  o'     —  o*  o  ° 

It  may  be  added  that  the  Laplace  transform  of  the  Green's  function. 


(2.11) 


00 


S'(r,0,sjr  )  =  Y2.      ^(r,e  +  2Bm,  sjr  )  , 


m=-oo 
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can  be  converted  into  a  Fourier  series  by  using  Poisson's  summation  formula. 
The  relevant  form  of  this  is 

00  .00  ,  Q 

(2.12)     T~      7(r,9  +  2iim,sir^)  =  ±^-75-  r_        F  (r,n,s}r^)e^"*'   j 


m=-oo 


(2n)-"/'  ^-00 


it  can  be  obtained  formally  by  calculating  the  Fourier  coefficients  of  the  left 

hand  side*.  Thus  e.g.  for  r  <  r  , 

^       —  o  ' 

00    r  l'(s)K  (sr^)K  (sr)  ,   .  „ 

(2,13)     G(r,0,8jr^)  -  £1   jy sr^)I^(sr)  -  ^ T-^ \  e^°« 

n=-oo  ^n 

which  can  be  written  as 

I*(s)K  (sr  )K  (sr)     00  l'(s)K  (sr  )K  (sr) 

(2.1U)  G(r,e,s}r  )  -  K  (sR)  -  ^ ?   °  ° 2  JZ_    -^ ?   °  " cos  n©  , 

°      °  K  (s)  ll^l      K  (s) 

o  n 

by  the  Addition  Theorem  for  Bessel  functions  of  imaginary  argximent;  K  (sR)  is  the 
Laplace  transform  of  the  elementary  solution  (l.lO).  Conversely,  (2.13)  and  (2.1U) 
co\iLd  have  been  obtained  first  smd,  using  Poisson's  summation  formula,  converted 
into  (2.11)  and  (2.10).  This  method  is  legitimate  but  less  transparent  than  that 
based  on  (1.12)  (which  is  actually  equivalent  to  (2.12)). 

The  integral  representation  (2.10)  of  a  typical  term  of  the  series 
(2.11)  is  similar  to  the  type  of  solution  obtained  in  the  discussion  of  the  pro- 
pagation of  electromagnetic  waves  in  a  stratified  atmosphere.   (In  fact,  the 
problem  of  finding  7  is  equivalent  to  such  a  case,  as  can  be  seen  by  making  the 
substitution  r  "  log  z) .  In  the  theory  of  propagation,  it  is  found  that  the  dis- 
turbance in  the  shadow  can  be  obtained  by  evaluating  the  integrals  in  question  by 
the  method  of  residues  in  temis  of  a  series  of  "propagation  modes"  (see  for  in- 

'  Sufficient  conditions  for  the  validity  of  Poisson's  formula,  in  a  different  fora, 
but  equivalent  to  (2.12),  are  given  by  E.G.  Titchraarsh,  Fourier  Integrals,  (Oxford 
1937)p.  60. 
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stance  Kerr'^-').  The  same  device  can  be  applied  to  (2.10).  For  the  integrand  is, 
for  fixed  s  /  0,  a  single-valued  function  of  v,  considered  as  a  complex  variable. 
It  can  be  shown,  for  s  real  and  positive,  that  the  contour  can  be  moved  across  the 
upper  or  lower  half-plane  of  v  according  as  ©  <  Oj  hence  (2«10)  can  be  replaced  by 
the  sxm.   of  the  residues  of  its  integrand  at  the  poles  in  the  upper  or  lower  half- 
plane  respectively.  These  are  given  by  the  zeroes  of  K  (s)  considered  as  a  function 
of  v»  It  is  shown  in  Appendix  A  that  for  s  real  and  positive  these  zeroes  are  pure 
imaginary;  hence  they  may  be  denoted  by 

(2.15)  V  -  itij     (J  =  1,2,...) 

'/  \ 
in  the  upper  half-plane.  As  K  (s)  is  even  in  v,  those  in  the  lower  half-plane 

jure  then  v  =  -ip...  Now 

J 

'l  (sr)K'(s)  -  l'(s)K  (sr)^  K  (sr  ) 

'  V      V  V     V     f   V    o 


is  also  an  even  function  of  j  j  thus  the  series  of  residues  equivalent  to  (2.10) 

is  formally 

l!  (s)K.   (sr  )K.   (sr) 
_  CO    in.    i|J..   o  i^. 

(2.16)     F(r,e,sir^)  =  2h  F"   J i- J exp(-R,|©|) 

[an  in  Jn  -  ^^ 

and  this  is  valid  for  both  r  <  r  and  r  >  r  .  Each  terra  of  (2.16)  satisfies  the 

—  o      —  o 

Laplace  transform  of  the  wave  equation,  and  is  the  Laplace  transform  of  a  '•  pro- 
pagation mode"  • 

3«  Approximation  Valid  for  Large  s 

When  s  is  large,  (2.16)  can  be  evaluated  approximately.  Since  the  be- 
havior of  F  for  large  s  determines  the  "initial"  behavior  of  F,  these  approxima- 
tions are  of  obvious  interest  in  the  theory  of  pvilse  diffraction.  By  hypothesis , 
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V  is  the  Laplace  transform  of  a  function  which  vanishes  for  t  <  0,  and  so  is 
regular  in  some  half-plane  Re(s)  >  c  >  0«  Hence  it  is  sufficient  to  take  s  as 
real  and  positive  in  the  first  instance.  Then  various  approximations  can  be 
used,  which  are  derived  in  Appendix  A. 

In  the  first  place,  if  s  — ♦  oo  and  the  order  j  is  fixed,  then 

1/3 

(3.1)  ti  ^  s  +  as  ^^ 

where  a.  is  the  j*^  zero  of  f  (a),   f(a)  being  defined  by 

This  can  also  be  written  in  terms  of  the  Mry  function*-  ■•, 

-°°  1/2 


/'"^  1/2 

cosC^Y^  +  Yz)dY     =     -YJT    \/A'^^^'^^' 


as 


(3.I4)  f  (a)      -     -^  Ai(-2^/^  a)        . 


Again,  if  from  now  on  |j..  is  understood  to  be  given  by  (3.1)  then 

J 

2^f\   f(a  )       ^ 

(3,5)  K   (s)  ^    YfY       exp(-^  nti^)        , 

J         3  s 


(3.6) 


2^/^n  a.f(a.) 


#-K.  (s)     ^  ^^2 J_   exp(-i-n^i.) 


An  approximation  for  K.   (sr)  valid  for  r  >  1  is  obtained  in  Appendix  A,  but 
when  it  is  substituted  in  (2.16)  the  result  is  too  complicated  for  the  calcu- 
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lation of  the  inverse  Laplace  transform.  One  is  therefore  forced  to  use  a 
simpler  fonnula,  valid  for  r  >  1  as  s  -)  co  which  is 


(3.7)   K^^  (sr)^ 

%i 


1/2 


n^  1„     1/2  ^a/2   -111    1/3   -11 


where  cos"  (l/r)  is  defined  to  be  between  0  and  -r  n.  For  r  «  1,  (3»5)  must  be  used. 
Also,  it  follows  from  (2.9)  and  the  definition  of  [i.   that 

l!   (s)  K.   (s)  -  i 
ijx.     1M-.        s 


whence,  by  (3«5) 
(3.8) 


J 


When  substituting  these  approximations  in  (2.16),  three  cases  must  be 
distinguished,  namely  r  >  1,  r  >  Ij  r  «  1,  r  >  1  or  r  >  1,  r  ■  Ij  and  finally 
r  =  r  =  1.  In  the  first  case,  a  typical  terra  of  (2.16)  becomes 


(3.9) 


U(r^  -  l)^^(r^-  1)V^ 


i-^  g /y   exp(-sf-a  5s  '■^) 
a,  liAo,;  Is'  *^ 


I'^^A' 


jL^  j 


where 
(3.10) 


lol      -11      -11 
|e|  -  cos    p  -  cos    —  , 


(3.11) 


X    '     5.   (r^.  1)^/2  .  (r^-  1)^2  ^ 


For  r  ■  1,  r  >  1,  one  has 
'  o    * 


(3.12) 


3"^^'  ^      exp(.sr.a^6sV3) 


2(r^-  1)1/^      V^V^ 


-12- 


where  6  .and  'X   must  rtow  be  evaluated  with  r  "  1»  The  typical  term  for  r  ■  1, 
r  >  1  is  obtained  from  (3 •IS)  by  replacing  r  by  r.  Finally,  if  r  ■  r  ■  1  then 

(3.13)  r  -  6  -  |el  (r  »  r^-  1) 

and  a  typical  term  is 

(3.1U)  -^^  exp  J-  |el(s  +  a^  s^/^) 

The  fractional  errors  of  these  approximations  are  all,  for  fixed  r  and  r  ,  of 
order  s~  '  • 

For  large  j  one  obtains  by  substituting  the  asymptotic  formulae  for  the 
Bessel  functions  in  (3*2)  that 

Henc%  as  j  -^  CO, 

(3.16)  a.~|(3n)2/3(j-^)^/'  . 

Thus  the  series  (3*9),  (3*12)  and  (3»ll4)  converge  only  if  6  >  0.  In  the  analogous 
problem  of  propagation  in  a  stratified  atmosphere,  which  has  already  been  referred 
to,  one  can  show  that  in  the  shadow  the  asymptotic  behavior  of  the  svun  of  the  series 
corresponding  to  that  for  F  can  be  inferred  from  the  behavior  of  the  individual  terms. 
Outside  the  shadow  this  is  not  the  case,  although  the  series  converges.  It  seems 
reasonable  to  assume  that  the  same  result  is  valid  in  the  problem  treated  here. 


U*  Approximate  Formulae  for  the  Green's  Function 

The  evaluation  of  the  inverse  Laplace  Transform  of  (3«9)>  (3»12)  and 
(3«lii)>  with  j  =  1,  involves  the  calculation  of  the  integral 
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C+ICD 


(U.1) 


Uj^(a^5,T)  =  A.  /     s"^  exp(sT  -  a^6  s^/^)ds 


C-lOO 


where 

(U.2)  T  =  t  -  "tr  , 

k  =  I",  ^  or  ^,  and  c  >  0,  This  integral  will  also  be  required  later  for  other 
values  of  k.  If  T  <  0,  the  contour  can  be  deformed  into  a  large  semi-circle  in 
Re(s)  >  0.  The  integral  taken  along  this  contour  vanishes  in  the  limit,  by 
Jordan's  Lemma,  Hence  U  vanishes  for  T  <  0. 

For  T  >  0,  U,  can  be  estimated  by  the  method  of  steepest  descents.  The 


function  sT  -  a, 6  '   of  s  has  a  saddle  point  when 

s  =  s^  =  (a^5/3T)3/2 
and  it  is  easily  shown  that 


/      2(a  6)^/^     a/2„5/2       2 
sT  -  (L,5  s^'^  "    -  -,A    ^  JA     +  ^ rrr-    (s-s„)  +  ... 


"1^  =      ■   -  pTT^   *  ^ 


W 


Kence  on  the  path  of  steepest  descent  through  the  saddle  point 

with  real  ^,  and  so 

(U.3)        U,(a,5,T)^  i^  (a,5)'3-6WA  T(6k-5)A  „(,)exp(-S^) 

where 


(i;.li) 


2(a^5)^/^ 

'1-  pyr^iT?    ' 


and  the  factor  H(T)  has  been  inserted  because  Uj^  "  0  for  T  <  0.  Using  (U.3),  one 
finds  the  following  approximations  for  F  from  (3»9),  (3.12)  and  (3.15):  • 
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1^/2      (r2-l)(r2-l)T6  V    L^^'^l^J 


8n- 

o 


(U.6)  F(l,e,tir  )  ^        g^       lA     W^ 

°  U(r>  l)-^/^  T^/^       a,f(a,) 

0  11 


(U.7)  F(i,e,tji) 


3^>^nl/2|e|^AH(T)        exp(-4i) 


l?7n ^ 


where,  in  (Ii.7)  T  =  t  -    |e|   and  6  =    |©|,  in  accordance  with  (3«lii)« 

If  3k  is  an  integer  then  U     can  be  evaluated  explicitly.     Take  first 

2  13 

k  =  -T  and  put  sT  =  ■!?■  sC^  in  (U«l)»     Then 

32/3,-1/3         />--P(W6) 

c/ooexp(-iTi/6)  *^  J 

Since  Re(s, )   <  0  for  r  "  <   |arg  s, |   <  ^  ">   the  contour  can  be  deformed  into  the 
imaginary  axis,  and  so 

-00 


ds^       o 


Thus  according  to  the  definition  (3 •3)  of  the  Airy  function, 

(U.8)  U2/3(a^6,T)   =  3^/^  T-^/^  H(T)Ai  L6(3T)-^M  -  \\\)         H(T)K^/3(C^) 

where  C-,  is  again  given  by  (li.Ii). 

The  value  of  U,  when  3k  is  an  integer  can  be  deduced  by  repeated  dif- 
ferentiation or  integration  with  respect  to  the  parameter  cu,6  from  (i;»8)  since 
clearly 
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U.9) 


dUj^(a^5,T) 


k^     -  -  ^3k-l)/3(^''^^ 


Hence  one  can  use,  instead  of  (U.5)  and  (U.7),  the  formulae 


/-\ 


(U.io) 


38/3  j^(^) 


Ai  } 


cu^6(3T) 


-1/3 


'^^''^'*^o)-^,(,2.,^lA(^..,)iA,V3         ;[,(^)]^ 


(U.ll) 


95-'-/^  H(T) 


K^(V 


U  n(r^-  l)lA(r^-  1)^^  ^V^   aj/^  [f(a^)]  = 


F(l,9,t;l)^-3^/W)  i  ^'  L5(3T)-1/3| 


-  H(T) 


V3^^]L^     ^  |Q|V3(gl) 


'^ 


3/2 1^,1/2  .1/2 


e 


3l7^T 


nl/3 


The  fractional  error  of  these  approximations  is  of  order  T  '  .  Finally,  the  cor- 
responding approximations  for  the  Green's  function  are  obtained  by  substituting  the 
appropriate  form  of  F  in  (1.13). 


5»  Interpretation  in  Terms  of  Geometrical  Optics 

To  interpret  these  results  it  is  best  to  begin  with  the  typical  constituent 
piilse  F»  This,  it  will  be  recalled,  is  the  field  due  to  an  incident  pvilse  defined 
on  a  Riemann  surface  whose  sheets  are  (2Tn-l)Ti  <©<  (2m  +  l)n  due  to  a  singularity 


which  has  the  character  of  an  instantaneous  source  at  time  t  =  0,  at  r  =  r  ,  8  =  0. 

'        o' 

If  the  cylinder  were  absent,  the  fronts  of  this  pulse  would  be  the  circles 
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(5.1)  t  =  R  »  (r^  +  r^  -  2rr  cos  o)^'^ 

at  least  for  0  <  t  <  r  .  The  fronts  of  the  secondary  disturbance  can  be  deduced  from 

—  0 

Fermat's  principle,  since  the  arguments  by  which  the  validity  of  this  principle  is 

[21 
established  for  pulse  fronts '-  -"  can  be  extended  to  the  case  of  propagation  of  a 

Riemann  surface.  There  is  a  family  of  reflected  fronts  E,  and  the  two  families  of 

diffracted  fronts  D  ,  D  •  The  fronts  E  are  the  normal  trajectories  of  the  rays 

obtained  by  reflecting  the  incident  rays  -  which  are  simply  the  straight  lines 

through  r=r,0=0-at  their  intersections  with  r  =  1  nearest  to  the  source. 
^      o' 

These  rays  issue  from  the  arc  |©|  <  cos"  —  ,  and  the  fronts  R  remain  in  the  sheet 

0 

ra  =  0,  For  0  <  t  <  r  -  1,  there  is  no  secondary  disturbance  at  allj  for 
r-l<t<(r-l)  '    ,  only  a  reflected  pulse  is  present}  but  for  *  >  (r  -  l)  '  , 
diffraction  into  the  shadow  takes  place.  On  the  Riemann  surface  the  shadow,  which 
is  bounded  by  the  prolongation  of  the  incident  rays  which  touch  r  =  1,  is  given  by 
the  inequalities 

(5»2)         9  >  cos"  -  +  cos"  —-  ,   ©  <  -  cos"  -  -  cos"  --  , 

r       r  '  r       r 

o  o 

or,  by  (3.10) 

(5«3)  5  >  0  . 


To  determine  the  arrival  time  of  the  diffracted  pulse,  and  the  diffracted 

I   "  J 

fronts  D  ,  D  ,  one  may  begin  by  considering  the  points  on  r  =  1.  One   front,  D  , 


-11  /  2  \l/2 

propagates  m  the  sense  of  increasing  0,  starting  at  0  =  cos   —  when  t  =  ( r  -1)  '    , 

It  propagates  onto  the  sheets  with  m  >  0.  The  other  front,  D  ,  starts  at 

-11         /  2  \l/2 
9  -  -cos   —  when  t  =  ( r  -11  '  and  encircles  the  cylinder  in  the  direction  of  de- 

0 

creasing  0,  propagating  into  the  sheets  with  ra  <  0.  Thus  by  (3.11)  the  equation  of 
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both  fronts  may  be  written  as 
($.U)  t  =  X 

so  that  t;  simply  denotes  the  arrival  time  of  the  diffracted  front  (only  one  dif- 
fracted front  at  most  arrives  at  a  point  on  the  Riemann  surface) *  The  rays  normal 
to  D  and  D   are  in  general  tangent  to  r  =  Ij  for  instance,  the  ray  normal  to  D 


through  (r,©)  touches  r  =  1  at  ©  -  cos   —  ,  and  is  of  length  (r  -  l)  '  •  Hence, 

o 


—  ,  and  is  of  length  (r  -  l)  '  , 

O  I        It 

again  by  (3«ll),  (5.ii)  is  also  the  joint  equation  of  D  and  D   at  all  points  of 
the  shadow.  Thus  T  is  the  time  counted  from  the  arrival  of  the  diffracted  front, 
and  the  factor  H(T)  in  all  the  approximations  obtained  in  Section  U  implies  that 
the  disturbance  starts  when  the  diffracted  front  arrives^  this  agrees,  as  it  shoixld, 
with  the  general  theory  of  diffracted  pulses. 
Now  G  is  given  by  (1.13), 

CO 

(5.5) 


G(r,e,tjr  )  -   r"   F(r,e  +  2mn,tjr  ) 


m=-oo  % 

Since  Q  has  period  2n  with  respect  to  0,  it  is  sufficient  to  take  0  <  |o|  <  n^  Let 

R        Inl        -11        -11 

Sq   =   |e|  -  COS    -  -  COS    —   , 

o 

(5.6)  6   =  2mn  +  e  -  cos   -  -  cos   -=-  ,      (ra  >  O)  , 

o 

5   =  _  2mn  -  9  -  cos   -  -  cos  ^     ,  (m  <  O)  , 

m  r       r 

o 

and 

(5.7)  1^.  "  5   -^  (r^-  1)^^^  *   (^o-  1^^'^^   '    T^  =  *  -  '^m  • 

mm  0  ra        m 

Consider  first  the  case  of  a  point  (r,e)  in  the  shadow,  so  that  5  >  0.  Then  the 

term  "  m"  of  (5.5)  is  zero  for  t  <  X   ,     Now  the  diffracted  fronts  in  the  physical 

—   m 

plane  are  obtained  by  •'  projecting"  D  and  d"  from  the  sheets  of  the  Riemann  sur- 
face onto  the  sheet  ra  =  D.  Diffraction  begins  with  t  =  (r  -  1)  '     at  the  two  points 


-  17  a  - 
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-1/  I   II 

r  =  1,  0  =  +  cos  (1/r  },  and  D  ,  D  ,  may  be  described  as  'coming  from'  these 

points  respectively.  As  t  increases,  the  end-points  of  D  and  D   travel  steadily 

round  the  cylinder  in  opposite  directions.  Thus  the  diffracted  fronts  meet  at  r  =  1, 

0  =  n  when 


t  =  n  +  (r^-  1)V2  _  ^^3-1  i 

o  r   ' 

o 

and  subsequently  penetrate  into  the  region  in  which  the  secondary  distiirbance  has 

already  been  established.  These  fronts,  and  the  reflected  ones,  are  shown  in 

Fig.  1. 

By  symmetry,  it  is  sufficient  to  consider  the  case  ©  >  0.  Then  the  first 

front  to  arrive  is  obviously  D  .  It  arrives  at  time  T  ,  and  G  is  then  given  by 

F(r,e,tjr  ).  Subsequently,  D  passes  (r,©)  when  t  =  "t  +  2mn  =  't  ,   (m  =  1,2,...). 

Each  time  that  it  does  so,  a  further  contribution  F(r,©  +  2mn,tjr  )  must  be  added 

to  G.  Similarly,  D   passes  (r,©)  when  t  =  X    ,   with  ra  =  -l,-2,,..  and  the  terra 

m 

F(r,©  -  2inn,tjr  )  must  be  added  to  G  on  each  occasion.  When  ©  <  0,  the  parts  played 

by  D  and  D   are  reversed j  on  ©  =  n,  both  arrive  simultaneously,  and  the  terms  with 

m  and  -m-1  in  (5*5)  are  equal  to  each  other. 

It  remains  to  discuss  the  case  of  a  point  (r,©)  in  the  reflection  region 

(5  <  0).   he  evaluation  of  (2.10)  by  residues  is  then  invalid  for  the  term  m  =  0 
0 

of  (5.5) •  The  other  tenns  still  represent  diffracted  pulses.  To  deal  with  the 
terra  m  »=  0,  one  can  apply  the  method  of  steepest  descent  to  (2.10).  It  is  found 

that  on  r  =  1,  F(r,©,t}r  )  is  simply  equal  to  the  doubled  elementary  solution  (l.lO) 

-1/2 
to  order  (t-R)  '  ,  which  is  to  be  expected.  One  may  therefore  identify  F(r,©,t5r  ) 

with  the  sum  of  the  incident  and  the  reflected  pulse  when  6  <  0.  The  calculations 

•^         o 

are  rather  involved,  and  do  not  furnish  any  new  information;  as  they  are  however  of 
some  analytical  interest,  they  are  given  in  outline  in  Appendix  B. 
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6.     Numerical  Res\ilts  for  the  Case  r  =  r  ■  1 

_^— — — — — ^-— ^— — — ^  o 

It  is  clear  that  the  approximations  obtained,  which  are  valid  when  T  is 
small,  can  only  be  applied  to  the  calculation  of  G  when  (5.$)  reduces  to  the  term 
m  «  0,  or,  in  the  neighborhood  of  e  «  n,  to  the  teims  m  =  0,  m  »  -1.  One  can  of 
course  also  estimate  the  additional  contributions  made  by  the  other  terms,  pro- 
vided that  T  is  small, 
ra 

Some  calculations  were  made  for  the  case  r  ■  r  ■  1,  which  is  probably 
typical,  using  the  value 

(6.1)  -^    a^/^  -  0.279   . 

Now  it  is  apparent,  from  (U»7)  and  (U.12i,  that  F  and  all  its  derivatives  with 
respect  to  T  vanish  for  T  «  0.  (In  fact,  F  has  an  essential  singularity  at  the 
diffracted  front,  and  this  is  the  reason  why  the  diffracted  pulse  cannot  be  ex- 
panded in  ascending  powers  of  T).  This  may  be  contrasted  with  the  behavior  of 

-1/2 

the  reflected  pulse  at  its  fronts  where  G  becomes  infinite  like  T  '     if  T  is 

the  time  counted  from  the  arrival  of  the  reflected  front.  However,  as  T  increases, 
F(l,©,t}l)  very  soon  increases;  and  as  6  --^0,  the  initial  rise  becomes  very  rapid. 
This  is  illustrated  by  the  curves  plotted  in  Fig.  2,  which  show  G  as  a  function 
of  T  for  ©«=n,  e-in  and  6  -  1^  n.  To  bring  these  curves  together  into  one 
figure,  a  different  scale  had  to  be  emplojred  for  the  ordinate  in  each  case.  Thus 
the  c\irves  are  the  graphs  of  10  G(l,  n,  t,  l),  0.5  G(l,  ^,  n,  t|  l)  and  of 
0.1  G(l,  r-,  n,  tj  1).  '^he  abscissa  is  T  =  t  -  ©. 
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FifTure  2 
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The  values  of  G  at  6  =  n  may  be  compared  with  those  of  the  elementary 
solution  G  given  by  (l.lO),  which  would  be  obtained  if  the  cylinder  were  absent, 
when  the  times  elapsed  since  the  arrival  of  the  respective  fronts  are  equal.  The 
front  of  the  elementary  solution  arrives  when  t  =  2,  and  with  T  ■  t  -  2  one  has 

(6.2)  G^  =  <)T*(T*+  h)[  H(T*)   . 


The  following  table  shows  the  values  of  G  -  calculated  by  means  of  (U»12)  -  and  of 

G  ,  for  0  <  T  <  O.U  and  T  -  T*", 
o'      —   — 

TABLE  1  (e  -  n) 


T  =  T 

0 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

O.UO 

G 

0 

0.03 

0.15 

0.27 

0.37 

0.U7 

0.53 

0.58 

0.63 

G 

CO 

2,22 

1.56 

1.27 

1.09 

0.97 

0.88 

0.81 

0.75 

Thus  the  cylinder  becomes  ineffective  as  a  ■  shield"  from  the  source  for 
pulses  whose  length  is  equal  to  about  ha3_C  the  radixis. 

For  small  9,  one  may  compare  G  with  the  solution  obtained  when  the  cylinder 

is  replaced  by  the  rigid  plane  x  =  1.  This  is  simply  2  G  ,  and  one  may  compare  the 

value  of  G  at  e  with  that  of  2  G  evaluated  at  a  distance  9   from  the  source.  This 

0 

gives  the  following  table  at  0  »  t-  n: 


TABLE  2 

(-^ 

n) 

T 

0 

0,05 

0.10 

0.15 

0.20 

0,25 

0.30 

0.35 

O.iiO 

G 

0 

5.50 

U.30 

3.65 

3.29 

3.01 

2.8U 

2.72 

2.61 

2G 
o 

00 

7.02 

Uo90 

3.9U 

3.37 

2.97 

2.67 

2.U5 

2.26 

Diffraction  of  Plane  Pulses 

The  diffracted  field  corresponding  to  any  incident  pulse  can  be  obtained 
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by  means  of  the  Green's  function.  As  an  example,  the  diffraction  of  a  plane 
piolse  will  be  considered.  Let  the  incident  pulse  be  given  by 


(7.1) 


p  ■  p  (t  +  x)  «  p  (t  +  r  cos  O) 


where  p  denotes  the  excess  pressure.  It  is  sufficient  to  take 


(7.2) 


p  (t  +  x)  »  6(t  +  x) 


where  6  denotes  the  Dirsic  6-f\inction.  For  if  P(r,e,t)  denotes  the  corresponding 
diffracted  field,  then  the  diffracted  field  is  given  in  general  by 


(7.3) 


p(r,e,t) 


P(r,©,t-t')p  (t')dt'  . 


CO 


Now  (7.2)  is,  formally,  the  inverse  Laplace  transform  of  exp(-sr  cos  9). 
But  the  Laplace  transform  of  the  elementary  solution  (l.lO)  is 


(7.U) 


G^(r,©,sjr^) 


K  (sR),   R  =  (r^+  r^-  2rr  cos  ^)^'^   , 

o     '  0      0* 


and  so,  for  large  r  , 


G(r,e,sir^)  =  (^   ^  e-^^  /l  +  0(4; 


°(i) 


2sr 


^-sro  ^sr  cos  eC  ^  ^^^.O  ^ 


It  follows  that 


(7.5) 

and  hence  obviously  also 


lim 

r  -^00 
0  ' 


'2='-  \'^'  sr  - 

°  '    e  o  G  (r,©,s) 


-sr  cos  © 


n 


(7.6) 


'2t  3 


1/2 


sr 


lim  {[•^  e  o  g'(r,e,s) 


r  -^  00 

0 


=  P(r,e,s)  . 


By  (1.13)  and  (2.16)  the  diffracted  field  P  is  therefore  given  by 
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00 

(7.7)  P(r,e,t)  =  21   Q(r,«  +  2nin,t) 

m=-cr) 

where  , 

I.  (s)K.  ,     . 
_  00    iti  '   iM..(sr) 

(7.8)  Q(r,e,s)  =  2nr"  jt^ A exp(-vL  |e|) 

j=l   I •5—  K.  (s) 


^=^J 


When  |s|  is  large,  Q  can  be  approximated  ty  means  of  the  formulae  (3.1)  to  (3.8) j 
alternatively,  the  approximations  (3.9)  and  (3 •IS)  can  be  substituted  for  F  before 
making  r  ->  00  in  (7.6).  It  is  found  in  this  way  that 


o 


(7.9)    Q(r,e,s)  ~  ^3/2  ^i/;^^;_  ^^1/1.    ^  r  ^  ^s^l/a   «P(-=6  -  a^Ss^S, 


when  r  >  1,  and 

(7.10)  Q(l,©,s)^  -i-       ^-^-^  exp  J-(|e|  -  I  n)(s  +  a^s^/^)  I 

where  6  is  now  obtained  from  (3.10)  by  letting  r  —^  00, 

(7.11)  5  -  |0|  -  in  -  cos'^  i 

and 

(7.12)  G-=   limCt+r)  -  5+  (r^-  l)"'"/^  . 


o 
r  -00 

0 


For  r  =  1 


(7.13)  5  =  cr  -  |e|  -  I  n  (r  -  1) 

The  series  (7.9)  and  (7.10)  converge  only  when  6  >  0»  This  again  corresponds  to 
a  diffracted  pulse,  for  the  shadow  boundaries  are  now  y  =  +  1,  'X  <  0  (they  are 
in  the  sheet  |©|  <  n  of  the  Riemann  siirface  on  which  Q  is  defined). 

When  6  >  0,  then  Q  vanishes  for  t<(j'}t"Cr'is  the  equation  of  the 
two  diffracted  fronts  on  the  Riemann  surface  corresponding  to  the  incident  front 
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t  =  -X  »  -r  cos  ©.  Thus  the  contribution  of  the  term  Q(r,S,t)  of  (7«7)  to  the 
Faltung  integral  (7«3)  is 

(7.11;)  /  p^(t  -  t')Q(r,e,t')dt' 


/ 


provided  that  p  (t)  -  0  for  t  <  0.  Since  the  approximations  for  Q(r,e,t)  which 
can  be  deduced  from  (7.9)  and  (7»10)  are  valid  only  when  t  -  C"  is  small,  (7»1U) 
can  be  applied  to  "  short"  pulses  only;  in  particular,  P  itself  can  be  considered 
as  giving  the  diffracted  field  of  a  short  "  square"  pulse  of  length  L,  if  it  is 
multiplied  by  L. 

8o  Evaluation  of  the  Diffracted  Field 

The  inverse  Laplace  transforms  of  (7»9)  and  (7«10)  are,  in  the  notation 
of  Section  U, 


(8.1) 


(8.2) 


Q(l,0,t)  ^  -i.   VAi^        , 
7^     a^f(a^) 


where  now 

(8.3)                    T  -  t  ~    Cr  , 

so  that  T  is  again  the  time  counted  from  the  arrival  of  the  diffracted  pulse.     Ac- 
cording to  (U.3),  (8.1)  and  (8.2)  give 

QfiV2   rr(^\  exp(-CJ 

(8.U                        Q(r,9,t)  ^  -^ '>         l/l  -TTTr ^ 

^3 A  fi3A  „rT^  exp(-4J 
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where  again  ,  /« 

2    (°1«^  ^ 


Also,  by  (U.IO)  and  (U.8), 


(8.7) 


VV'^>  -  ^  Al"|c^6(3T)-V3J  =  ^^  Ai|;^6(3T)-V3J 


so  that  one  may  replace  (8«5)  by 

(8.8)        Q(l,e,t)  ^  J^^^^  ^  V3^^1^  • 


53/2  jj(^)  a^/ 


As  in  the  case  of  the  Green's  functions,  the  series  (7»7)  reduces  to  a 
finite  sum.  Each  term  is  a  contribution  which  only  differs  from  zero  after  the 
passage  of  one  of  the  diffracted  fronts  across  (r,6).  In  particular,  if  one  takes 
^  n  +  cos~  —  <  ©  <  n,  then  the  first  pulse  to  arrive  is  Q(r,6,t)«  Its  front  pro- 
pagates in  the  positive  sense  round  the  cylinder,  and  at  its  m-th  passage  the  term 
Q(r,e  +  2(m+l)n,t)  must  be  added  to  P.  The  other  set  of  pulses  associated  with  a 
front  propagating  in  the  ncigative  sense  round  the  cylinder,  is  given  by  the  terms 
of  (7.7)  with  ra  <  0.  At  ©  =  n,  the  terms  with  m  and  -m-l(m  >  O)  axe  equal  to  each 
other.  The  total  excess  pressure  P  can  only  be  calculated  by  means  of  (8.1i)  and 
(8.5)  or  (8.8)  immediately  after  the  first  arrival  of  the  diffracted  front}  one  can 
however  say  that  in  general  the  additional  contributions  due  to  the  terms  with  m  >  1 
or  m  <  -2  are  negligibly  small. 
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When  the  values  of  a^,  f(a^)  are  substituted  in  (8.1;), (8. 5),  these 
equations  become 

5,  exp(-^) 

(8.9)  Q(r,e,t)  ^    3.5U  5 y^    y^  H(T)  , 

4^^  exp(-^;, 

(8.10)  Q(l,e,t)  ^    8.16  -i = 2^  H(T)  , 


with 


iV^ 


(8.11)  ^     -  0.28  -^ 


T 
Some  curves  of  P  against  T  computed  from  these  formulae  are  shown  in 

Pig.  3  and  U*  Figure  3  gives  (P,T)  curves  for  r  =  1  and  6  =  0  -  -^  n  equal  to 
30  ,  hy  f   60  and  90  j  Fig.  U  gives  the  (P,T)  curves  on  ©  =  n  for  the  same  values 
of  6,  which  correspond  to  r  =  2,  l.hDi,  1.155  and  1.0.  Both  sets  of  curves  show 
that  for  T  >  0.05  P  is  roughly  independent  of  5  and  rj  they  also  illustrate  how 
the  diffraction  effect  lessens  as  one  moves  away  from  the  point  r  =  1,  ©  =  n  either 
along  the  cylinder  or  along  the  line  ©=n.  Atr=l,  e=n  the  infinitely  short 
incident  pulse  is  converted  into  a  pulse  consisting  of  a  roughly  linear  pressure 
rise  followed  by  a  very  slow  decrease,  so  that  diffraction  appears  to  have  an  in- 
tegrating effect. 
9.  The  Diffraction  of  Short  Harmonic  Waves 

The  diffraction  of  harmonic  waves  can  be  treated  by  the  simple  device  of 
putting  s  =  ixo  in  the  Laplace  transform,  where  \   =  2n/oo  is  the  ratio  of  the  wave- 
length of  the  incident  waves  to  the  radius  of  the  cylinder,  and  -multiplying  by 
e"""  .  For  instance,  if  the  incident  waves  are  plane,  so  that 

(9.1)  pjt)  -  e^*  , 

then  (7.3)  gives 


40- 


30- 


20- 


10- 


-  25  a  - 


8=  30« 


3  =  45* 


Figure  3 


.06 


-  25  b  - 


8- 


4- 


r  =  2 


Vs  1.414 


r=l.l55 


r  =  l 


.02 


.04 


.06 


Figure  L 
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P(r,e,t-t')  e^^  dt' 


1  I 

but  as  PCr,©,t-t  )  =  0  for  t  -  t  <  CT"  this  is  equivalent  to 

(9.2)  p(r,G,t)  -  e"^*/  P(r,©,t')e"^*  dt'  =  e"^*  P(r,e,i^)  . 
Hence  uQ- 

.  .    00     _ 

(9.3)  p(r,e,t)  =  e    ^Z    Q(r,e  +  2ran,ioo)  . 

in=-oo 

Each  term  ^(r,©  +  2inn,  ico)  of  this  series  is  a  sum  of  certain  modes,  obtained  by- 
putting  s  ■  ico  in  (.78).  As  co  -^oo,  the  asymptotic  approximations  can  be  used. 
For  instance,  when  r  =  1,  a  typical  mode  is  asymptotically 

1/2 


(9.1i)    -^   ^ 


V^  a.  f(a.) 


.[l©  *  2mnl  -  IJ  [i(.  .  >  .V3)  .  3^  ,^^^/^ 


>   • 


This  is  a  simple  wave,  propagated  with  phase  velocity  (1  +  ^  a.o)"  '  )"~  and  subject 

to  an  exponential  attenioation  which  reduces  the  space  of  the  amplitwie  by  l/e  in 

-1/2  -1  -1/3 
the  distance  3    a.  oj  '  •  The  attenuation  increases  both  as  co  and  the  order 

increases;  also,  if  co  is  large  then  the  terms  with  m  >  0,  ra  <  -1  in  (9.3)  are 

negligible.  Because  of  the  phenomenon  mentioned  at  the  end  of  Section  3,  only  the 

modes  J  =  1  can  be  used  to  estimate  p  for  lairge  oo.  The  resulting  expressions  are 

p  -  1.837  exp  i  ^  [t  -(1+0.118x2/^) (©  -  J)]  -  1.295x"^/^(e-J) 


+   1.837  exp 


r^  [^  -(1.0.118x2/3)  (|1  .  e)J.  i.295x-^/\^  -  «) 


-27- 


Near  e  =  n  these  two  terms  are  of  comparable  amplitudes,  and  as  |0  -  n|  increases 
IpI  has  maxima  and  minima  approximately-  ^  X  apart j  but  as  |©  -  n|  increases  fur- 
ther, one  of  the  terras  dominates  the  other  one  and  |p|  increases  exponentially. 
When  X  is  small  enough,  one  may  define  an  effective  shadow  boundary  as  the  value 
of  e  -  ^  or  as  -w-  -  e  at  which  |p|  equals  a  given  small  number  e,  and  assxirae  that 
only  one  of  the  terms  in  (9»5)  needs  to  be  taken  into  accoxmt.  The  distance  of 
this  effective  shadow  boundary  from  the  geometric  shadow  boundary  is 


(Q  A\  Q   "  -  >■     1   1.837 

(9.6;  e  -  ?;■  =  r^o^  log  — ;^- 


.V3 


that  is  to  say  this  angle  varies  as  the  cube  root  of  the  wavelength. 
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Appendlx  A 
Bessel  Functions  of  Iinaginary  Order  and  Argument 

I, In  this  appendix,  an  outline  of  the  derivation  of  the  approximations 
quoted  in  Section  3  will  be  given.  The  results  are  applied  to  Laplace  transforms, 
which  sire  necessarily  regular  in  sone  half -plane  Re(s)  >  c  >  Oj  hence  it  is  suf- 
ficient to  take  s  to  be  real  and  positive* 

The  differential  eqioation  satisfied  by  K.  (s)  is 


2 

(A.i)  k:;(s)  .  i  k:^(s)  ^(^-  i)k.^(s)  =  o . 


Prora  this  and  the  corresponding  equation  obtained  by  replacing  [i  by  its  complex 
conjugate  p.  it  follows,  unless  |a.  is  real,  that 

(A.2)    /'  Vs)K^^*(s)  f  =  -^   k^x*(-)<^.(^^  -  V^-)V^^ 

Now  K.  (s)  has  the  integral  representation'-  '  ^*   -" 

00 

(A. 3)  K.  (s)  »  /   e""     ^  cos  [ly  dy  • 

Hence  K.  ^(.(s)  is  the  complex  conjugate  of  K.  (s),  ajid  the  left  hand  side  of  (A, 2) 

is  essentially  positive.  It  follows  that  neither  K  (s)  nor  K.  (s)  can  vanish 

for  complex  ^,  as  this  would  imply  the  vanishing  of  K.  *(s)  and  K.  »(s)  respectively 

and  render  (A.2)  contradictory, 

I 

Thus  the  zeroes  of  K.  (s),  which  are  required,  occur  for  real  |j.  only, 

when  s  is  real.  Let  p.  be  fixed,  and  s  varied.  When  s  is  large  then 

V     f   \      ,   f  ^  \l/2  -s   ^  v'    f    \  c  "  ^l/2  -s  <  0  . 
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Hence  as  s  is  decreased  from  +00,  K.  (s)  is  positive  and  increases  at  first.  Let 

s  be  the  greatest  s  such  that  K.  (s)  »  0.  Then  K.  (s  )  is  positive  and  a  maximum, 

so  that  K.  (s  )  =  0.  But  by  (A.l) 

2    2 

K.  (s  )   =  ^  K,  (s  ) 

iM.  o        s      Im-  o 

and  so  necessarily  \\i\  >  s   *  This  inequality  remains  true  if  s  is  fixed  and  11  is 
variedj  hence  the  positive  zeroes  of  K.  (s),  |j.  ,  |j.p,  ..•,  all  satisfy  the  inequality 

(A.U)  ^ i  -  ^   (j  =  1»  2,...)  . 

n.    When  3  is  large  then  by  (A.U)  the  \i.  are  also  large,  and  K.     (s)  can 
be  approximated  by  Langer's  method'^-'   (a  convenient  sunmary  is  given  by  Lighthlll '-  -* ) . 

Let 

(A.5)  u(z)      =     K^i^(tiz)      , 

when  z  is  real  and  positive.  Then 


(A.6)  u  (z)  +  iu  (z)  +  /(  ^  -  l)u(z)  =  0  . 

z 

Hence  for  large  p.,  and  z  <  1, 

(A.7)  u(z)  ^  ^^  S)1A  ^  V3^^^  *  '  '^-1/3^^^]  ' 

where  A,  B  are  independent  of  z,  and 

(A.8)  K 


/  (  4  -  1)^^^  dZ  -  IX  cosh"^(l/z)  -  n(l  -  z^)^/^  . 


z 

To  determine  A  and  B,  let  z  -♦O.  Then  on  the  one  hand 


V^'  ■  riECiSCT  [w^' -  V^^)]  ~  jfe  |(i^)'(|^-)-*''-(-i^)!c|^-)^ 


whence 
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(A,9)  u(z)  .-^  (-J-)  '   exp(-^  iiti)sin(M.  log  ^  +  ^)  • 

On  the  other  hand,  X. '^  '^  log(2/ez)  as  z  •>  0,  and  so  (A, 7)  gives 

ii(z)  .^  (-)  '     J  A  cos(p.  log  ^  -  r  -  j)  "^  B  cos(ii.  log  ^  -  r  +  |) 

By  comparing  this  with  (A, 9)  one  finds,  after  a  short  calculation,  that 

n  exp(-i  n^i)^^/^      p 
(A.IO)  u(z)  ^  ^1/2  ^l/2(.^  _  ^^)1A  |V3^^^  "  -^-1/3^^^  !>  • 

Now  let  z  ■  s/m-j  so  that  ^  is  defined  by 

(A.ll)  X     '     V-  cosh"-"-  (n/s)  -  {\?  -   s^)-"-/^  J 

then 

n  exp(-i  ^^t)^^/2  r 
(A.12)  K.^(s)  ^  ^iy^^^2_  ^^)1A  ^1/3^^^   *   ^-1/3^^) 

The  fractional  error  of  this  approximation  is  0(|j."  '    ),  and  since  p.  >  s  by  hypothesis, 

—2/3 
this  can  be  replaced  by  0( s     ) • 

III.  Let  s  -^  CO  in  such  a  manner  that  ^  rernains  finite.  Then  "C/s   -^0, 
and  by  (A.ll)  this  implies  that  |i/s  ->1.  When  (m./s)-1  is  small  then 

where  the  dots  stand  for  a  series  in  ascending  powers  of  (ti/s)-l.  Hence  (ti/s)-l 

2/3 
is  a  series  in  ascending  powers  of  (^/s)  '  .  Let  therefore 


(A.IJ4)  ti  =  s  +  as-"-/^  +  oCs"-"-/^) 

Then,  by  (A,13) 

(A.15)  ^  =  3-  a^/^  +  0(s"^/^) 


and  also 
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^2  -  s^     =     2as^/^  .  0(s2/^)      . 


When  these  approximations  ai^  substituted  in  (A, 12)   and  in  the  expression  by  differ- 
entiating (a, 12)  with  respect  to  s,  one  obtains 


(A.16) 

(A.17) 

where 

(A.ie) 


v^^^ 


2^/2 
3s 


/2  C  " 

^  exp(-|  nn)f(a)  Jl  +  o(s"^/^) 


,1/2. 


V^=^    "   -  7^  ^^"^  "t^)f'('x)  ii  +  o(s~2/^) 


3s' 


f(a)     =     a 


1/2 


J         ,  23A  a3/2   \                 /23/2  ^3/2 
V3l 3 j    ^  V3( J 


>  • 


Hence  the  zeroes  of  K.  (s)  are  given  by 


(A.19) 


ji.  -  s  +  a.s-^/3  ^  o(s"-^/3) 


where  ou.,  a^f  •••  are  the  positive  zeroes  of  f  (a),  arranged  in  ascending  order. 
This  is  true  as  s  — >  oo  for  fixed  jj  the  constants  implied  by  the  O's  depend  on  j, 
and  presumably  tend  to  infinity  with  j. 

IV.  It  is  now  possible  to  derive  approximations  for  the  functions  reqiiired 
in  the  evaluation  of  (2.16).  In  the  first  place,  one  may  differentiate  both  sides 
of  (A, 17)  with  respect  to  a   and  then  put  \i   «  \i.»     This  gives,  because  of  (A.lii), 


(A.20) 


=1/3^ 


[|r  <.(=']...j  -  -  ^  -p(-i  ""j'^'^-j)  [i  *  °(="'/')]    • 


Now  f(a)  can  be  related  to  the  Airy  function  by  the  equation 


(A.21) 


f(a)  =  -^  Ai(-2^/3  a) 


and  so 
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(A,22)  f     (a)      -     2ar(o)      . 

Hence  (A. 20)   can  be  replaced  by 

3/2 


(^•23)  [i:  </^) 


.../  ^  v^v  J^  *  °^^"/^^ 


Again,  K.      (sr)    is  given  by  (A. 12)  with 


(A.2U) 


, -1  n^       ,   2         2  2.1/2 


For  fixed  r  >  1  and  s  -*  oo,  (x.  <  sr  xiltimatelyj  hence  it  is  convenient  to  replace 


L.12)   by  an  alternative  formula.     Now  for  p.    <  sr. 


(A.25) 


.3in/2, 


2  2       2.1/2 


e'-'^y\,.,         ix,.     =     (s  r  -  p..)  '     -  cos      —      , 


r^lii 


sr 


and  so  (A. 12)  cein  be  written  as 


(A.26) 


K...   (sr) 


expC-rr  nii  )7\, 


1/2 


iTT^j 


..^^Vsr;  ^    ^^2^2_  ^.^-lyi^.  .^/3V'.j 


(n,J  il  +  0(s'^/^)i  . 


This  is  equivalent  to  (A. 12)  and  either  form  is  valid  for  all  r  >  1.     But  as 
(A.26)   is  too  complicated  for  the  evaluation  of  the  inverse  Laplace  transform 
of  a  typical  mode,   one  must  approxiinate  further,  taking  r  >  1.     Then  for  large  s 


u. 


r.   2     ,.1/2      -111      1/3   -1  1  ^  ^/  -1/3n 
'(r  -  1)  '  -  cos  -S  -   a.s  '  cos   -  +  OU    )  , 


and  as  this  makes  7]^.     large,  one  can  put 


Hence  (A,26)  can  be  replaced  by  the  cruder  approximation 


(A.27)  K.  (sr)^ 


ili 


■J 


(2s) 


n^/^       _^  fl       r  2  ,.1/2    -lll^    ,, 

i/2(^^_.,^)i/it  ^p  <r2  ""^f  ^L^^  -^^    -^°^   ^^'^j^  °°^ 


1/3   -1  1 

ri 
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Appendlx  B 
The  Reflected  Pulse 

When  the  quantity  6  defined  by  (3»10)  is  negative,  the  series  ex- 
pansions obtained  by  approximating  the  fornnila  (2.16)  for  F  do  not  converge, 
and  it  may  be  assumed  that  the  evaluation  of  (2,10)  by  means  of  residues  - 
which  yields  (2,16)  -  is  not  s\ii table  for  large  s.  Some  other  method  of  evaluat- 
ing (2,10)  must  then  be  adopted.  The  simplest,  and  mosf  satisfactory,  case  is 
r  =  1.  Then  (2,1j)  reduces  to 

/>°°   K  (sr  )   ,  ^ 
(B.l)  F(l,©,s}r  )  =  -  /     -n~^  e^       dv  , 

/co    ^V^) 

Now  when  s  is  large,  real  and  positive,  then  K  (sr  ) A  (s)  is  roughly  inde- 
pendent of  V  when  |v/s|  is  small j  hence  the  principal  contribution  to  the  in- 
tegral (B,l)  arises  from  large  values  of  |vl,  and  one  may  attempt  to  evaluate 
it  by  substituting  appropriate  asymptotic  expressions  for  K  (sr  )  and  K  (s). 
These  can  be  derived  by  applying  Jeffreys'  method  (the  WKB  method)  to  the  dif- 
ferential eq-oation  of  K  ,  or  from  (A, 26),  Dropping  the  suffix  and  using  the 
asymptotic  formula  for  K^  /^,  one  can  write  this  as 

(5-  n)  '  -L 

V^  ^  7T-2 — ZTTn  ^p(-  ?  "  t^  -n. )  . 

^  (s  r     -  |i  )   ' 

Now  let  |j.  "  ivj  then 

U  -     (s'r^.  vy/2  .  iv  cos-^(.  1^)     =     (sV.  v2)l/2_  ,  ,^^^'1  ^  .  1  ivn  , 

and  so 

(1  ^)l/2  p  -| 

(B,2)  K^(sr)  ^  —j-g 2  yi^     exp  Jv  sinh""^  ^  -(sV+  v^)'^^^\ 
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-1  V 

This  is  Valid  for  complex  v  if  the  imaginary  part  of  sinh   —  is  suitably 

restricted;  if  in  particular  v  is  pure  imaginary  (and  |v|  <  sr)  then 

-i  sinh"  —  =  sin"  (-  ~)  must  be  between  -  -rr  n   and  +  ^  ii«  ^-Iso 
sr  sr  2        2 

(B.3)        sK^Cs)  ^  .(I  n)V2(32,  ,2>>lA  ^  T  3^^-!  ^  .(32,  ^2^^1A|  ^ 
When  (B,2)  and  (B,3)  are  substituted  in  (B.l)  one  obtains 
(B.li)  F(l,e,s}r^)  ^  /   e^"^  g(v)dv 

(y-  CD 

where 


(B.5)        h(v)  =  ive  +  V  sinh"^  iF"  -  '^  sinh"^  ^  +(s^+  v^)-^/^-(sV+  v^)"^/^  , 

o 

(B.6)  g(v)  =  (s^.  v2)-lA(sy.  v^)-^^  . 

One  can  now  evaluate  (B.U)  approximately  by  the  method  of  steepest  descents. 
The  function  h(v)  has  a  saddle  point  at  v  »  v  ,  where 

(B.7)  ie.  +  sinh   —  -  sinh   _£  -  0  . 

sx*  s 

The  solutions  of  this  equation  can  be  obtained  geometrically,  as  follows.  Let 
A  be  the  point  (1,©)  (where  F  is  to  be  calculated)  and  B  be  the  point  (r  ,0) 
(where  the  soiurce  is  situated) .  If 

(B,8)  a  =  zlOBA  and  p  =  n  -  ZlOAB,  then  6  =  p  -  a,  sin  p  =  r  sin  a  . 

Hence  (B«7)  can  be  solved  by  putting 

,  V  ,   V 

(B.9)  sinh   —  "     -iB,  sinh   — ^  =  ia  , 

3      "^ '        sr 

o 

and  because  of  the  restrictions  |p  |  <  »■,  |a|  <  *■  n,  this  is  the  only  admissible 

solution*  Also,  one  must  have  |g|  <  cos   —  that  is  to  say  6  <  0» 

0 


-35- 

At  the  saddle  point, 

(B.IO)  h(v)  =  s(cos  P  -  r  cos  a)  =  -  sR 

where 

(B.ll)  R  =  (r^  +  1  -  2r  cos  g)""-/^ 

is  the  distance  ABj  also 

(Ba2)   h"(v)  -  -  3r  cos\  cos  p  >    ^^V  "  /    ^^^  \  ^^^   ,\l/2    * 

o  sir  cos  a  cos  p)  ' 

Since  h  (v  )  is  real  and  negative,  the  line  of  steepest  descent  through  the  saddle 

point  is  parallel  to  the  real  axis(it  is  assumed  that  s  is  real  and  positive),  and 

hence 


s^r  cos  a  cos  Py 

The  inverse  Laplace  transform  of  this  is 

(B.U)       F(i,e,t}r^)  ^  J  5^~rrr  I     "^*  -  ^^  • 

Now  the  elementary  solution  G  ,  defined  by  (l.lC),  is  approximately  given  by 


1/2   .     ^ 

t  -  R) 


(B.15)  Go(l,e,t  jr^)  ^  |^(t\R)T  '  H( 

at  (l,e),  and  so  (B.lU)  can  be  written  as 

(B.16)  F(l,e,t;r^)  ^  2G^(l,Q,tjr^)   . 

But  this  is  eqiiivalent  to  the  well-known  result  that  the  pressure  is  doubled  at 
the  pulse  front,  on  a  reflector.  One  may  therefore  conclude  that  F(l,©,tjr  )  re- 
presents the  sum  of  the  incident  and  reflected  pulse. 

If  r  >  1,  and  6  <  0,  one  cannot  proceed  immediately  in  this  manner*  For 
it  is  now  necessary  to  separate  the  incident  and  reflected  pulses.  It  is  not  pos- 
sible to  do  this  at  once  with  the  expression  (2.10)  for  F,  which  is,  for  r  <  r  , 

X   -,       ^     /'°°  Vsr)K'(s)  -  l'(s)K  (sr)        . 
(B.17)   F(r,G,s}r^)  -  /   -^^ ^ ; ^ i^—  K  (sr  )e^''®  dv  , 


-00  ^v^^^ 
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while  r  and  r  must  be  interchanged  when  r  >  r  .  For  the  resulting  integrals 
o  —  0 

would  be  divergent.  But  one  may  write  (B.l?)  as  the  sum  of  two  cosine  transforms, 

p^                                                      /.°°i;(s)K  (sr)K  (sr  ) 
(B.18)  F(r,0,sjr  )  =  2  /   I  (sr)K  (sr^)cos  vedv  -  2  /   ^-5 cos  v©dv  . 

°    /o  Jo  \^'^ 

The  first  term  is  obviously  the  elementary  solution  defined  by  (l.lli)  on  the 
Riemann  surface,  when  the  cylinder  is  absent.  Since  this  is  the  solution  generat- 
ed by  an  instantaneous  source  at  r  =  r  ,  9  =  0,  it  follows  from  the  uniqueness 

theorem  that  it  will  be  identical  with  G  for  0  <  t  <  r  j  when  t  >  r^,  there  will 

o  0—0 

be  an  additional  distiurbance  whose  front  is  t  ■  r     +  r.     Thus  it  will  be  equal  to 

G     for  small  values  of  t  -  R  at  any  point  of  the  reflection  region,  and  the  second 
o 

term  in  (B»18)  will  be  the  additional  disturbance  due  to  the  presence  of  the  cy- 
linder.    This  then  is  itself  a  solution  of  the  wave  equation,  and  reduces  to 

F-G^G     onr  =  l:itisa  diverging  pulse  which  vanishes  for  t  <  R  on  the 
00  ' 

cylinder  r  =  1,  and  it  must  therefore  have  the  reflected  front  as  its  pulse  front. 
It  is  therefore  the  reflected  pulse.     As  a  matter  of  fact,  if  the  Bessel  functions 
are  replaced  by  their  asyngstotic  approximations,  two  integrsuLs  are  obtained  each 
of  which  can  be  estimated  by  the  saddle-point  method,  and  it  is  found  that  the  sura 
of  the  contributions  from  the  two  saddle  points  is  the  leading  tem  of  the  express- 
ion of  the  reflected  pulse  in  half-integral  powers  of  the  time  elapsed  since  the 
arrival  of  the  refle  cted  front. 
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different  from  the  one  presented  here. 
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Errata 

Page  12:  Replace  the  paragraph  after  Equation  (3.16)  by  the  following: 

The  question  which  now  arises  is  to  what  extent  both  the  series 
(2.16)  and  the  approximations  (3.9)  ^  (3.12),  and  (3.II4)  can  be 
used.  This  has  not  been  discussed  in  the  present  problem,  but 
will  be  treated  in  the  paper  on  propagation  in  a  stratified  at- 
mosphere which  has  already  been  referred  to.  The  results  obtain- 
ed there  suggest  that  as  s  — >oo,  F  can  be  replaced  by  the  approxi- 
mate expression  {3  •?) ,   (3.12)  or  (3.lU)j  for  the  first  mode  (j=l) 
of  the  series  (2.16)  provided  that  5  >  0.  Consequently,  attention 
will  from  now  on  be  confined  to  this  term  only. 

Page  13:  Second  line  of  second  paragraph: 

read   sT  -  cu,6s  ^^     instead  of  sT  -  a,  6  '  . 

Page  18:  Sixth  line  of  second  paragraph: 

read     F(r,0+2mJT,t;r  )         instead  of  F(t,@  -  2mn,t;r  ) 

Page   29 J      Third  line   in  first  paragraph: 

read  K.    (s   )  <  0  instead  of  K.    (s   )   =  C 

111     0'  i|j.^    0' 

fourth  line  of  first  paragraph: 
read 


2         2 
^0  -^ 

instead  of 

2         2 

^     -  % 

s 
0 

Page  30:  At  the  end  of  the  paragraph  before  'III'  add  the  following: 

The  order  terms  given  here  and  in  III  may  differ  from  the 
correct  ones.  However,  for  the  applications  in  this  paper  these 
differences  are  inconsequential. 
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